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The regime of small isospin chemical potential in QCD is investigated. Using the phase quenched 
partition function in the e-regime an expression for the chiral condensate is given, which is studied 
in the temperature isospin chemical potential plane. Lines of constant values of the condensate are 
shown and it is estimated how the critical temperature varies as a function of the isospin chemical 
potential. Finally, the dependency of the fermion sign problem on the chemical potential and 
temperature is examined. 
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I. INTRODUCTION 



During the past years the phase diagram of QCD 
at finite baryon density has been a subject of great 
interest. Eventhough, a large number of phases have 
been suggested [11, Q only now direct confirmation by 
first principle calculations has emerge. One reason is 
that the phase of the fermion is complex, which makes 
standard Monte-Carlo simulations possible only for 
small values of the chemical potential [1, 0, H, H, 0,11, [1] . 
The phenomenology of heavy ion collisions and neutron 
stars has also been contributin g t o the confirmation 
of the phase diagrams of QCD [10]. In both neutron 
stars and relativistic ion collisions the isospin density is 
different from zero. Therefore, it is of great interest to 
study the influence of isospin chemical potentials on the 
phase diagram of QCD. The main goal of this paper is 
thus to examine the phase diagram of QCD at nonzero 
temperature and isospin chemical potentials. 
The starting point will be the work of Gasser and 



Leutwyler, who propounded that in the e-regime [11 1 
the low energy partition function, which is dominated 
by the zero modes of the Goldstone bosons, reduces to 
a group integral uniquely determined by the pattern 
of chiral symmetry breaking [l^, [H, [Tj]- Since the 
lightest degrees of freedom, the pions, have zero baryon 
charge, the low energy effective partition function at 
zero temperature is not affected by the baryon chemical 
potential. In addition, due to the complex fermion 
determinant it is not possible directly to simulate lattice 
QCD at nonzero baryon chemical potential by proba- 
bilistic methods. By instead considering the absolute 
value of the fermion determinant, one gets a situation 
which is doable. For an even number of flavors this 
corresponds to the product of a fermion determinant and 
its complex conjugate. The conjugate flavors correspond 
to ordinary fermionic flavors with the opposite sign of 
the chemical potential 15]. A theory with one flavor and 
a conjugate flavor is therefore identical to a theory with 
two flavors at nonzero isospin chemical potential ^iso- 
For further discussions of isospin chemical potential see 
[TgI . [TtI [TsI . IToj . Since the pions have nonzero isospin 
charge, the low energy effective partition function de- 



pends on the chemical potential. Thus, by looking at the 
isospin chemical potential instead of a baryon chemical 
potential, a situation is reached, where numerical calcu- 
lations and nontrivial analytic computations are possible. 

The organization of this paper is as follows. In section 
// the expected features of the QCD phase diagrams for 
nonzero temperature, baryon and isospin chemical po- 
tentials are reviewed from general arguments. QCD is 
considered in the e-regime for which the volume of Eu- 
clidean space-time is taken such that chiral perturbation 
theory is valid and that the Goldstone modes associated 
with chiral symmetry are the dominant degrees of free- 
dom [lO]. Moreover, the Goldstone field is treated as 
constant. This has the advantage of allowing for exact, 
analytic calculations. In section /// the effective theory 
is used to derive an expression for the chiral condensate 
in the case of two quark flavors Nf = 2. The expression 
is presented in the (/ii^o, T)-plane and from the lines of 
constant values of the chiral condensate it is estimated 
how the critical temperature drops as a function of the 
isospin chemical potential. In section IV the fermion 
sign problem is investigated by plotting the expectation 
value of the phase factor in the (/x, T)-plane. Finally, in 
section V concluding remarks are given. 



II. QCD AT NONZERO CHEMICAL 
POTENTIAL 



The QCD partition function at temperature 1//3 and 
chemical potential is given by 



Zqcd = X! ^ 
fc 

where the sum is over all states. This partition function 
can be rewritten as the Euclidean QCD partition function 
for Nf quark flavors 

Nf 

Z''^■{m■^i)^{\{dci{D + ^Jino + mf)), (2) 

/=1 

where the average is over the Euclidean Yang-Mills ac- 
tion and the Dirac operator is D = 7^i(9^ -I- i^^t), with 
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7;i the Euclidean 7-matrices, and an SU{Nc) valued 
gauge potential. The quark masses are denoted hy ruf 
and the chemical potential for each flavor is denoted by 
fif. Below, focus will mainly be on QCD with two quark 
flavors and nonzero baryon number and isospin chemical 
potential, which in this case are defined as 





fJ-B = 


T 


f^iso — 




Deconfined 




<'i''P>=0 


Confined 




<ii"P>+0 





(3) 
(4) 









Deconfined 




<^kp>=0 


Confined 




<^>¥) 






/ Confined 




<JC>+0 



111^/3 



m. 12 



FIG. 1: A schematic phase diagram of QCD in the tempera- 
ture baryon chemical potential plane (left) and the temperature 
isospin chemical potential plane (right). The tUtt denotes the 
pion mass and mjv denotes the nucleon mass. 

Before proceeding to the expression for the chiral 
condensate some general remarks on the QCD phase 
diagrams are given. For a further discussion of the QCD 
phases see |21, 22, 23]. Figure [T] shows the simplified 
phase diagram of QCD in the temperature baryon 
chemical potential plane and in the temperature isospin 
chemical potential plane. At temperature T ^ 170 
MeV QCD undergoes a phase transition from a chirally 
broken phase to a chirally symmetric phase. Since lattice 
QCD requires a real action, information of the structure 
of the phase diagram along the chemical potential axis 
cannot be achieved reliably using that regularization. 
Therefore, most of the knowledge of QCD at nonzero 
baryon density is based on models and gene ral p roperties 
of the phases of QCD. For a discussion see [24| . 
In the case of isospin chemical potential an analysis 
of the chiral Lagrangian to one-loop order in the low 
temperature range shows that the pions Bose con- 
dense for fiiso > m7r/2, where m^r is the physical pion 
mass [2^, [2^. Contrary to the partion function for 
7^ 0, the partition function for fiiso 7^ can be 
simulated by Monte-Carlo methods. When using an 
isospin chemical potential instead of a baryon chemical 
potential the fermion determinant is real and therefore 
traditional numerical lattice methods can be applied [2^ . 

Measurements of lattice QCD are necessarily per- 
formed at finite volume. In this paper the standard 
choice in lattice regularization of QCD is used. This 
choice considers a torus, where the extension in the 
Euclidean time direction determines the temperature 



L4 = 1/T and the quantities Li,L2,L3 specify a 
three dimensional box. By considering a torus with 
Li = L2 = = L4 = L the volume of the Euclidean 
space-time is thus defined by = L''. 
Moreover, the effective theory is considered in the pres- 
ence of an isospin chemical potential. With the usual 
pattern of spontaneous chiral symmetry breaking for 
two light fiavors, the theory is described by a Lie-group 
valued field U{x) £ SU{2). As mentioned, focus will be 
on the e-regime, also known as the microscopic domain 
of QCD. In this regime the chiral Lagrangian is treated 
as a perturbative expansion around the zero-momentum 
modes in a finite volume. 



III. EFFECTIVE THEORIES AT LOW 
ENERGIES 

In QCD the low-lying degrees of freedom are the pions, 
which according to Goldstone's Theorem result from the 
spontaneous breaking of chiral symmetry for Nj > 2. 
Since the mass of the pion is about 140 MeV, pions are 
much lighter than the lightest non-Goldstone particles, 
such as the nucleons, which has a mass of about A ^ 
1 GeV. Therefore, at sufficiently low energies the QCD 
partition function is well approximated by the partition 
function of an effective low-energy theory involving only 
pions. The QCD partition function in a finite Euclidian 
volume V4 = is dominated by the pions if 



(5) 



This statement follows by comparing the contribution 
of the pion, exp(— to^L), to that of a heavier par- 
ticle, exp(— AL). The interactions of the pions can 
be separated into zero-momentum modes and nonzero- 
momentum modes. It was realized by Gasser and 
Leutwyler that there exists a kinematic regime, where 
the fluctuations of the zero-momentum modes dominate 
the fluctuations of the nonzero- momentum modes [3l| . 
This regime is given by the condition 



1 



(6) 



where is the pion mass. Intuitively this means that 
the wavelength of the pion is much larger than the lin- 
ear extent of the box. Thus, the pion field does not 
vary appreciably over the size of the box, which result 
in small derivative terms. Therefore, it is possible to 
only consider the zero-momentum modes in the e-regime 
In this regime the Lie group can be factorized 
as U{x) = J7o cxp[iV2'/>(a;)/F^], where Uo is the zero- 
momentum part and (j){x) represents the fluctuation fields 
[3^ . It turns out that to leading order one keeps only the 
static modes in the path integral, while the fluctuation 
degrees of freedom decouple. This can be seen as fol- 
lows. When the isospin chemical potential couples to the 



3 



theory it gives rise to a covariant time derivative in the 
effective SU{2) Lagrangian [H, [H, [13] , 

doU{x) ^ VoU{x) = doU{x) - /i,,o[^3, U{x)], (7) 

where is the usual PauU matrix. The leading-order 
terms in the effective Lagrangian then read 

C = ^Tr[VoU{x)VQU\x)+d^iU{x)diU\x)] 

-^Tr[MU\x)+M^U{x)], (8) 

where M. = diag(TOM, is the quark mass matrix and 
E is the chiral condensate. By expanding the Lie group 

U{x) - [/o[l + iy/24){x)/F^ + •••], (9) 

the usual kinetic term for 4>{x) is produced. When 
power counting is used in the e-expansion it is assumed 
that m^r ^ — 0{e^), while (j){x) ^ 1/L = (e) and 



that a consistent power counting for the /Uiso-term is 
IJLiso = C'(e^) 0, EU, [13]. Indeed, when the covariant 
derivative Eq. ([7]) is expanded using Eq. ([9]), the leading 
contribution becomes 



Vot/(a;) = iV2/F^do(j>{x) - ^Ji,so[<J3. C/o] + • • •■ (10) 

In the chiral Lagrangian Eq.® the mixed terms 
dQ4>{x)[ay,, Uq] produce only boundary contributions and 
play no role here. Thus, to leading order in the e- 
expansion the fluctuation field 4'{x) gives rise only to the 
kinetic energy term [33| 

j d^x^Trd^cj,{x)d^^{x), (11) 

which decouples as in the theory with /Ziso = 0. 
Collecting the remaining terms it is seen that the leading 
contribution to the partition function in the e-regime is 
the zero-dimensional integral 



Z^/=2(X;^,,„) = / rf[/eiVF,VLoTr[C/,<T3][C/t,.3] + iSVTr(^tc;+^C;t)^ ^^2) 

JSU(2) 

where the O-suffix on the group element U G SU (2) has been dropped for convenience. 

Projection onto fixed gauge field topology v is done by a Fourier transform, and amounts to the simple modification 

Mm 

Z^f=\M;^i^.o)= [ d^7(detC/)'^e3^^"^•-^'■[^■'^^l[c/^-3l+iEVT.(A1t^/+M£/t)^ ^^3^ 

where the leading-order contribution to the e-regime depends only on the scaling variables 

m, = TO^EF, ALo = ^^loF^V. (14) 
Derivation of the integral Eq. (|13p gives the following expression [s^ 



2 /■-^ 2 2 

{mu.rhd-, fiiso) = 2e^'^"° dt t e~^''"°* /^(i 'mu)Iu{t rhd). 
Jo 

In sectors of fixed topological index ly the chiral condensate follows by differentiation of Eq. 



where for simplicity to = rhu = rhd- 

By summing over topology Eq. (|16p gives 



^ /o^dte-^AL*- E:::-oo^.(^<)(w(mt)+/._i(TOt)) 

J^dt e-2ALot' (/i(2TOt) +/„i(2TOt)) 



(15) 



ii^i^ju ^ -drh^ogZ^' [m; Hiso) ^ — T]— 77^7^ — 7-— , (16) 

2 dt t e-^t^i^ot^I^imt)^ 



(17) 
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For simplicity the vacuum angle 9 has been set equal to zero 3£ 
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FIG. 2; The chiral condensate is plotted as a function of the 
isospin chemical potential for m — 1. The upper line shows 
the full chiral condensate The lower line shows the case for 
1/ = 0. It is seen how the simplified case, v = Q, makes up a 
fine candidate for the full situation. 

In Fig. [5]Eq. PT)) is plotted as a function of ^liso- The 
upper line shows the case where a summation over v has 
been made, and the bottom line shows the case where 
V = Q. It is seen that the simplified case, v = does 
not diverge significantly from the full situation. Thus u 
is set equal to zero from now on and the expression of 
the chiral condensate is given by 



2j^dt e-2ALot' ^2 i^(^rnt)h{rnt) 



Io{mt) 



(18) 
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the massive particles becomes increasingly important. 
Thereby, the theory cannot be described by only con- 
sidering the zero momentum modes in the e-regime. 
By performing a three-loop analysis of the low temper- 
ature region of QCD, Gerber and Leutwyler have esti- 
mated the value of the critical temperature at which the 
chiral phase transition takes place This is done by 
plotting the temperature dependence of the chiral con- 
densate and following read off the temperature in the 
limit where (i^ip) = 0- Their analysis indicates that in 
the chiral limit the phase transition occurs around 

^ 170 MeV. 

This value is beyond the range of validity of the formulas, 
but since the order parameter falls rapidly at the upper 
end of the range, Gerber and Leutwyler consider it mean- 
ingful to make the estimate. 

In this paper the critical temperature is estimated as the 
maximum value of the differentiated condensate with re- 
spect to temperature for fixed values of /Ltiso- 
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FIG. 4: The critical temperature, Tc, is plotted as a function 
of the isospin chemical potential, fiiso, for {2)j,^F^)/Tq = 10. 
The figure shows how Tc slightly drops with increasing ^lao- 

In Fig. 0] the critical temperature is plotted as a 
function of /i^so < m-^ll. It is seen how Tc slightly drops 
for increasing isospin chemical potential. 



FIG. 3: The figures show the evaluation of the chiral con- 
densate in the {^iso,T)- and {^1so,T)-plane respectively. The 
value of the condensate decrease with darker colours. To per- 
form the plot one free parameter has to he chosen. In order to 
make the values on the axes dimensionless the free parameter 
has been picked to he {2ficF^) /Tq — 10, where /ic = It 
is emphasised that To is not the critical temperature at which 
the QCD phase transition takes place, but is simply the scale 
chosen for these particularly plots. 

In Fig. [3] lines of constant values of the chiral 
condensate (tpip) are plotted in the {fj,iso , T)- and in 
the (/if^o , T)-plane respectively. The increase in the 
imaginary plane is in agreement with [37| . 

The validity of the formulas fails at a temperature of 
T ^ 150 MeV. At this point the interactions among 



By Taylor expanding the critical temperature, 



Tcilliso = 0) 



= l-Ai 



Hi. 



(19) 



A calculation of the average value of the curvature for 
{2^lF^)/Ti^ = 10 gives A « 0.01. This value is relatively 
stable for {2ijl^F^) /T^ > 5 and is comparable to the value 
obtained in [27|] for the lowest quark mass. 



IV. THE FERMION SIGN PROBLEM 

As mentioned, the fermion determinant is complex 
at nonzero chemical potential. This prohibit standard 
Monte Carlo sampling of the path integral in what is 
known as the sign problem. Overviews with references 
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to various studies of this problem may be found in 

UMMM- 

By factoring the determinant of the Dirac operator into 
its absolute value and the phase factor, exp^icjj), the de- 
terminant can be rewritten as 

det(i:> + m + ^7o) = I det(L» + m + M7o)|e*'^. (20) 

The severity of the sign problem can be measured 
through the expectation value of the average phase fac- 
tor. A physical interpretation is obtained by defining 
the phase with respect to the phase quenched partition 
function 




1 2 3 4 5 



FIG. 5: The expectation value of the average phase factor, 
{e^'''')pq, is plotted in the {fj.,T)-plane for {2ixiF^)/T^ = 1. 
The value of {e?^"^) pq increase with lighter colours. Since the 
fermion sign problem is mild, when the phase factor is of the 
order one, the figure shows how the sign problem is getting 
more mild with increasing temperature. 

According to ^Ol] the sign problem is mild if (e^*"^)pq is 
of the order one. In Fig. [5] Eq. (|^ is plotted for = 
in the (/x, T)-plane. A plot after summation over v has 
been made too, but since no significantly difference be- 
tween the two plots were to be seen, only the case for 



= is shown. Thus, the v = Q contribution makes up 
a fine candidate for the visualization of the sign problem. 
The figure shows that the value of (e^*"^)pq increase with 
larger temperature and thus makes the sign problem 
mild. Figure[5]is in good agreement with ,41], where ran- 
dom matrix theory (RMT) is applied to obtain an ana- 
lytic expression for the average phase factor as a function 
of T, hb and m. 

V. CONCLUSIONS 

In this paper QCD has been considered in the e-regime 
of small isospin chemical potentials. Contrary to the 
regime of baryon chemical potentials, where the complex 
fermion determinant makes it impossible to directly sim- 
ulate lattice QCD by probabilistic methods, the isospin 
chemical potentials makes up a very accommodating sit- 
uation. 

The e-regime was specified such that chiral perturbation 
theory was valid and that the Goldstone modes associ- 
ated with chiral symmetry were the dominant degrees of 
freedom. Moreover, the Goldstone field was treated as 
constant. This had the advantage of allowing for exact 
analytic calculations. 

Further, the effective theory has been used to derive an 
expression for the chiral condensate in the case of two 
quark flavors Nf = 2, which has further been plotted in 
the {^iso,T)-\Aa.iie and in the (/il^^j, r)-plane. The lines 
of constant values of the chiral condensate showed a cur- 
vature in the real plane and an increase in the complex 
plane. 

Additionally, an estimate of the critical temperature as a 
function of the isospin chemical potential has been given. 
The measured small decrease of the critical temperature 
proved to be comparable to the values obtained in lattice 
QCD [27^. 

Finally, the fermion sign problem has been investigated 
by plotting the expectation value of the phase factor in 
the (/i, T)-plane. It was seen how the phase factor gets 
larger with increasing temperature and thus makes the 
fermion sign problem more mild. 
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